Non-Abelian spin singlet states of two-component Bose gases in artificial gauge fields 
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We study strongly correlated phases of a pseudo-spin- 1/2 Bose gas in an artificial gauge field 
using the exact diagonalization method. The atoms are confined in two dimensions and interact via 
a two-body contact potential. In Abelian gauge fields, pseudo-spin singlets are favored by pseudo- 
spin independent interactions. We find a series of incompressible phases at fillings v = 2fc/3. By 
comparison with the non- Abelian spin singlet (NASS) states, constructed as zero-energy eigenstates 
of a (fc + l)-body contact interaction, we classify the non-trivial topology of the states. An additional 
spin-orbit coupling is shown to switch between NASS-like states and spin-polarized phases from the 
Read-Rezayi series. 

PACS numbers: 73.43.-f,03.65.Vf,37.10.Vz 

Keywords: Gauge fields. Ultracold Bose gases. Correlated states. Exact diagonalization. 



Topological phases of matter attract a great deal of 
attention since the late 1980s, when the fractional quan- 
tum Hall (FQH) effect of electrons in two dimensions 
(2D) subjected to a strong perpendicular magnetic field 
was recognized as a paradigm of such phases (l[ . It was 
then pointed out that amongst these states there might 
be some which support non-Abelian excitations, where 
the exchange of two quasiparticles is described by some 
unitary matrix, as it shifts between different degenerate 
ground states Ql- Recently, this property has triggered 
the interest in topological phases, as quantum gates op- 
erating with non-Abelian anyons niight pave the way to 
fault-tolerant quantum computers [3|. For this purpose, 
the possibility of controlling or even designing quantum 
many-body states seems to be advantageous if not in- 
evitable. A very important tool to do so are cold atoms 
which can, in various ways, be manipulated through in- 
teractions with light Despite being electroneutral, 
the atoms may even couple to external gauge potentials, 
which, instead of being true electromagnetic vector po- 
tentials, have to be synthesized. Artificial magnetic fields 
have been achieved experimentally by rotating the atomic 
cloud [HI (see also @ for a review), though instabilities 
at high rotation frequencies have so far prevented from 
reaching the highly correlated regime. A promising alter- 
native pursues the generation of gauge fields by a laser 
dressing of the atoms [7|. This method has already been 
proven experimentally [s. It is readily generalized to 
non-Abelian gauge fields [9| , and has recently allowed for 
synthesizing a spin-orbit coupling [lo| . 

The strongly correlated states occurring in such sys- 
tems have been numerically investigated for the case of 
spinless or spin-polarized bosons , in which the incom- 
pressible phases are well described by the Read-Rezayi 
(RR) series of parafermionic states [1^, including the 
Laughlin wave function [l^, and the Moore-Read (MR) 
wave function Although also bosons with a pseudo- 
spin-1/2 degree of freedom have already been condensed 



14| , and their behavior under rotation has been studied 
both theoretically [3] and experimentally [l^ , few atten- 
tion has yet been paid to the strongly correlated regime 
in such systems. For the electronic FQH effect, however, 
a variety of states have been proposed for systems with 
spin, which despite the strong magnetic field might be 
in a singlet phase due to the small gyromagnetic ratio of 
electrons in many solids. Within a conformal field the- 
ory (CFT) framework, a generalization of the RR series 
to the so-called non-Abelian spin singlet (NASS) series 
[l3l has been formulated. Here the term "non-Abelian" 
refers to the nature of the excitations. For bosons, this 
series of topological states occurs at fillings v = 2fc/3, 
and includes the Halperin state [13] with Abelian excita- 
tions at fc = 1. Like the RR states, the NASS states are 
constructed as the zero-energy eigenstates of a repulsive 
[k + l)-body contact interaction. 

In this paper we show, using exact diagonalization, 
that also two-body contact (2BC) interaction leads to the 
formation of incompressible states precisely at the NASS 
filling factors in a 2D Bose gas exposed to an Abelian 
gauge field. For fc = 1. this is trivially the Halperin 
state, but also for k ~ 2 significant overlap with the cor- 
responding NASS state and a similar spectral structure 
indicate that the system can be brought into the NASS 
phase within this realistic setup. 

This finding is in close analogy to the spin-polarized 
case with RR-like states [ll| at filling v = k/2. However, 
as it is experimentally difficult to get into the regime of 
low filling, the two-component system benefits from sup- 
porting analogs of the RR states at higher filling. Fur- 
thermore, its possibility of tuning independently inter- 
spin and intra-spin interactions provides additional con- 
trol. For pseudo-spin-independent repulsive 2BC interac- 
tions, a mechanism is discussed which favors pseudo-spin 
singlet states. We show that, counterintuitively, states 
with zero interspin interaction energy are obtained by 
increasing the interspin interaction strength. Finally, ef- 
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fects stemming from a spin-orbit (SO) coupling can be 
integrated in the framework of two-component gases by 
generalizing to a non-Abelian gauge field. We investi- 
gate this scenario, and find a transition from the single- 
component RR states 1^ 2^ to the NASS states at SO 
coupling strengths corresponding to a Landau level (LL) 
degeneracy. 

We consider a two-component Bose gas in 2D, which on 
the single-particle level is described by the Hamiltonian 
^sp = (p ^ ^)^/(2"^)- The vector potential A is given 
by 



B(0,a;,0)l + g(f7^,cry,0). 



(1) 



The first term, acting equally on both pseudo-spin com- 
ponents, describes a magnetic field of strength B in the 
Landau gauge. The second term, with ai being the Pauli 
matrices, accounts for a Rashba SO coupling controlled 
by the coupling strength q, which we set to zero in the 
first part of this work. A detailed description how to im- 
plement this Hamiltonian for electroneutral atoms can 
be found in the supplementary material of Ref. [3] ■ As 
the boundary conditions we choose a periodic square of 
length a. This eliminates edge effects from our study, 
and, moreover, generates a non-trivial toroidal geometry, 
on which topological phases can be classified by the num- 
ber of degenerate ground states Solving the single- 
particle Hamiltonian yields, for q = 0, the usual LL 
structure [2l|, where every state |n,j, s) within a fully 



degenerate Landau level n is characterized by a momen- 
tum quantum number j and an additional two-fold de- 
generate spin quantum number s J,. Since an energy 
gap 2B separates the different LLs, we restrict ourselves 
to the lowest LL, n = 0. Fixing the number of fluxes 
per unit cell, iV$, restricts j to values running from 1 to 
N^, and the many-body Hamiltonian is, up to a constant 
amounting for the kinetic energy, given by 



H=En 



Si ^^232 %3S3 %4S4 I 



(2) 



where a,js annihilates a particle in state |0,j, s), and 
{j, s} denotes the set of quantum numbers ji, . . . , j4 and 
si, . . . , S4. The interaction conserves total momentum 
J = mod iV$, and the spin of the particles, so the 

coefficients read 



(3) 



(E) {j2,S2\gsiS25{zi - Z2) \j3,S3) |j4,S4) , 

with gsiS2 the spin-dependent contact interaction 
strength of two particles, and 6' a Kronecker delta mod- 
ulo N^. The Hamiltonian is diagonalized in blocks with 
a fixed particle number N, defining the filling factor 
1/ = iV/iV$, fixed Haldane momentum K [i^l, and a fixed 
spin polarization S = — N^, where (-^i) denotes 
the number of spin-up (spin-down) particles. 



First; we show that for spin-independent interactions, 
i.e. 5siS2 = 9^ the singlet state is energetically fa- 
vored. This configuration is a natural choice, as many 
experiments lJ| work with the \F = 1, uif = 0, f ) states 
of 23Na, or the \F = l,mj^ = -1) and \F = 2,mF = 1) 
states of ^^Rb with almost equal s-wave scattering 
lengths within and between the components. We start 
by contrasting the singlet to the fully polarized situation, 
5 = A^, in which we recover the results from Ref. [lH, 
with incompressible phases at v = fc/2, forming the RR 
series. The fully polarized zero-energy state with highest 
filling has = 1/2, and is the Laughlin state. In the 
opposite limit of a fully unpolarized system, zero-energy 
states occur up to ly = 2/3, where the 221-Halperin state 
is the exact unique ground state of the 2BC interaction. 
The analytic expression of this state on a disc reads 



i<j 



k<l 



i,k 



where we omit the irrelevant exponential term 
exp ( — X) I •^i P /4) ■ This wave function vanishes whenever 
two particles are at the same position, and it is symmet- 
ric under exchange of two spin-up or two spin-down par- 
ticles. However, it is antisymmetric under exchange of 
particles with different spin, and thus cannot be used for 
describing spinless bosons. For this case, also the terms 
in the last product had to be squared, which would result 
in a less dense Laughlin-like wave function. Generalizing 
this observation, we state that any spatial wave function 
which solves the fully polarized problem must also be a 
solution for \S\ < N, while the opposite is not true. Thus, 
for a spin-independent interaction, it strictly follows for 
the ground state energies E{Si) with 1511 < |5'2| that 
E{Si) < £^(5*2). Our numerical data, partially shown in 
the left side of Fig. [U supports this finding. In fact, 
in all cases we have investigated, the ground state of the 
S = 2 subspace occurs in the spectrum of S* = as a low- 
lying excitation, except for A^ = 6 and A^j, = 8, where 
we have a degeneracy, E{S = 0) = E{S = 2). In most 
cases, thus, the ground state is a singlet, which is un- 
derstood by noting that interactions of a ti-pair can be 
fully avoided by a single zero of the form Zi-^ — Zki in the 
wave function, while at least two zeros must be spent if 
a pair of particles with equal spin is to be anticorrelated 
in a similar way. 

Feshbach resonances allow to tune the spin- 
dependency of interactions (cf. Q). Increasing g-^^ 
while leaving g^^ = gn = g should favor spin polar- 
izations with less t4-"Pairs. However, this does not 
necessarily drive the system into the fully polarized 
configuration: As shown in Fig. [2] for A^ = 6 at i/ = I, 
the energy of even the fully unpolarized system saturates 
for large g^^. In this limit, the real ground state has 
5 = 2 and = 0.449 (in units of gB), close to the 5' = 
and 5 = 4 ground states with E = 0.466 and E = 0.475. 
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FIG. 1. (Color online) Left: Ground state energies, not show- 
ing possible degeneracies, at A*"* = 6 for different A'^ and 5* 
subspaces with gsis2 = 9- Even (odd) values of S correspond 
to even (odd) N . Right: A vs. u for N,s, — 6 and gsi32 = 9- 
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FIG. 2. (Color online) Interaction energy for A'^ = 6 and 
5 = 0,2 as a function of g-fi- Left: A'^* = 6. Right: A'^s = 5. 



All of them are well separated from the fully-polarized 
ground state with E = 0.581. In contrast to this, the 
energy of spin-unpolarizcd systems will never saturate 
at larger filling, as is shown in Fig. [5] for v = 6/5. This 
different behavior is understood by noting that for ly = 1 
the number of available zeros in the wave function is 
sufficient to completely suppress interactions between 
pairs of different spins, (E^i) 0, while it is not for 
larger ly. Jumps in the curve of (E-fi) as a function of 
g-fl show that the states with (Efi) = are reached 
by several abrupt re-organizations of the ground states. 
We note that once the state with (E^^) = is obtained, 
this state must be an eigenstate for arbitrary g^^. Thus, 
by abruptly switching off the interspin interaction, one 
could produce binary mixtures of highly entangled, 
non-interacting systems. 
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FIG. 3. (Color online) Energy spectra of singlet states for 
TV = 12 particles at filling v — 4/3. Black crosses (with the 
left energy scale) show the results for an SU(2)-symmetric 
2BC interaction, while red triangles (with the right energy 
scale) are the NASS spectrum for any three-body interaction 
with 5sis2S3 > 0. We do not show the three-fold center-of- 
mass degeneracy, but visualize further degeneracies by slight 
shifts in K. 



Coming back to spin-independent interactions, we ask 
now the question whether, in analogy to the RR states 
for spin polarized systems, incompressible states can be 
found amongst the spin-unpolarized ground states. Since 
incomprcssibility is connected to a discontinuity in the 
chemical potential, we define the particle-hole excitation 
gap as A{N) = N[E{N + 1) /{N + 1) + E{N - 1) /{N - 
1) - 2E{N)/N], which is plotted for iV* = 6 in Fig. [U 
(right). Upwards peaks correspond to downward cusps in 
the energy as a function of z^, and signal incompressible 
phases. They occur at fillings v = 2k/2> with fc G N. For 
these fillings, a series of singlet states being exact zero- 
energy ei gens tates of a (fc -I- l)-body contact interaction 
is known 



17| . For fc = 1, this is the 221-Halperin state, 
which has Abelian excitations, while the states for k > 
1 are predicted to have non-Abelian excitations, which 
established the name NASS series. 

We first construct these states for fc = 2 and 3 by di- 
agonalizing the three- and four-body contact interaction 
for up to = 12 particles: As predicted by CFT, these 
states are characterized by a (fc -I- l)(fc -I- 2)/2-fold topo- 
logical degeneracy on the torus. A subtlety occurs at 
fc = 3, where the CFT construction [13, being restricted 
to A^ = 2kp with p G N, is not able to predict zero-energy 
ground states for all A^ possible at v = 2. In fact, we find 
zero-energy states also for A^ = 8 and 10, but, in con- 
trast to the states with A^ = 6 and 12 which have the 
predicted ten-fold degeneracy, they are non-degenerate, 
and thus belong to a different topological phase. 

The overlaps of these states with the ground states of 
Eq. ([2]), given in Table U turn out to be significant: For 
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TABLE I. Overlaps between zero-energy eigenstates of (k + l)- 
body interaction with corresponding eigenstates of the 2BC 
interaction at energy i?2b (in units of gB) and Haldane mo- 
mentum K (in units of 2TTh/a). All states at fiUing f — 4/3 
have an additional three-fold center-of-mass degeneracy. 



k = 2, the states are exactly equal for = 4, and even 
for a relatively large-sized system with = 12 the over- 
laps are still around 0.8. Such a decrease of the overlap 
with large N also occurs in the case of the RR states in 
single-component systems with 2BC interaction Q, sug- 
gesting the use of ground state degeneracies to charac- 
terize the topological phase in the thermodynamic limit 
Topological degeneracies are usually lifted in finite 
systems, as also happens here. In the case of = 4/3 we 
may consider the fact that the six lowest eigenstates of H 
agree in their Haldane momenta K with the ones of the 
sixfold degenerate three-body eigenstates as a signal for 
belonging to the same topological phase. Both spectra 
are shown for = 12 in Fig. [3l We note that both the 
degeneracy splitting Aqs and the gap to the excited band 
Aex behave non-monotonic with N. The situation is dif- 
ferent ioT v — 2, where the candidates for the degenerate 
ground state manifold mix up with excited states. 

Our results therefore suggests that only the u = 4/3- 
NASS phase can be produced with 2BC interaction. For 
the experimental feasibility of this state to be high, it 
should also be robust with respect to deviations from the 
SU(2)-symmetric interaction. In the case of a (fc -I- 1)- 
body contact interaction at v = 2fc/3, the numerical val- 
ues of each interaction parameter do not influence the 
ground states. Though this is different for 2BC interac- 
tion at = 4/3, the overlaps with the NASS states re- 
main high within a wide range oi g-^^/g, and have a max- 
imum in the SU(2)-symmetric configuration. For N = 8, 
all overlaps arc above 0.8 if g\i/g S [0.8, 1.6]. 

Finally, we consider an additional SO coupling with 
q > in Eq. ([T]). This coupling mixes internal and 
external degrees of freedom on the single-particle level, 
such that the LLs become superpositions of the form 
\n,j,±) = at \n -IJA)+Pn where and 

are functions of q [19|. For general q, this lifts the previ- 
ous spin-degeneracy, as the single-particle energies now 

I. As has been 



read = 2Bn + 2q^ ± ^/W 



pointed out in Rcfs. [3, E^, the system then behaves 
basically like a single-component system, and the incom- 
pressible states are described by the RR series projected 
into the lowest LL of the SO coupled system. However, 
for q^ /B = 2n -I- 1, a degeneracy between E~ and 
occurs as the single-particle ground state switches from 
jn, 7, — ) to \n + l,j, — ). The Hamiltonian from Eq. ([2]) 
can then be adapted to the system with SO coupling by 
identifying the s quantum number with the two-fold de- 
generacy between n and n + \ , and by re-defining V{js} 
accordingly. Note that in this case 5* does not correspond 
to the spin polarization, but quantifies the population 
imbalance between [n, j, — ) and |nH- l, j, 

As observed in Ref. [20[ within a lowest LL approxima- 
tion, the projected RR states do not describe the system 
close to these degeneracy points. For q^ / B = 3, the sys- 
tem becomes gapless at both the Laughlin and the MR 
filling V ^ 1/2 and 1. However, we find a clear gap in 
the spectrum at v = 2/3 with iV = 6, and a small gap 
at 1/ = 4/3 with = 8. Projecting this = 2/3 state 
into the lowest LL of the system without SO coupling, 
and comparing with the 221-Halperin state, we find an 
overlap of only 0.19. It should be noted that, in contrast 
to the previous case of a pure spin degeneracy, interac- 
tions now can induce fiips within the two-fold degenerate 
manifold, and therefore S is not conserved. For a com- 
plete description of the state, we use generalizations of 
the Halperin state, obtained by diagonalizing the 2BC 
interaction at i/ = 2/3 in an Abelian field {q ~ 0) within 
different spin subspaces. Then the highest overlap of 0.60 
is obtained with the state having S = —6, which corre- 
sponds to all particles being in |0, j, — ). This is also the 
most populated LL, showing that, despite the degener- 
acy on the single-particle level, interactions favor the LL 
with lower n. The total fidelity for -6 < 5 < 6 is 0.82. 
For I' = 4/3 with A^ = 8, but now comparing the non- 
dcgcncratc ground state of the SO coupled system with 
the ground states of a three-body contact interaction at 
q = for different S*, we still obtain a total fidelity of 
0.52, summing from —8 < 5* < 8. 

We have investigated two-component Bose gases in ar- 
tificial gauge fields. The possibility of tuning the inter- 
action constants allows to control and manipulate the 
strongly correlated nature of the eigenstates. Our main 
finding is the formation of incompressible phases de- 
scribed by the NASS series within the experimentally 
feasible setup of particles interacting via two-body con- 
tact interaction. The internal degree of freedom permits 
also to realize systems with spin-orbit coupling. 
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Note added: After submission of this article, a preprint 
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by S. Furukawa and M. Ueda appeared presenting addi- 
tional numerical evidences for the NASS phase [23[. 
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